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ABSTRACT 

We demonstrate a new technique for determining the physical conditions of the broad line emitting 
gas in quasars, using near-infrared hydrogen emission lines. Unlike higher ionisation species, hydrogen 
is an efficient line emitter for a very wide range of photoionisation conditions, and the observed line 
ratios depend strongly on the density and photoionisation state of the gas present. A locally optimally 
emitting cloud model of the broad emission line region was compared to measured emission lines of four 
nearby (z « 0.2) quasars that have optical and NIR spectra of sufficient signal-to-noise to measure their 
Paschen lines. The model provides a good fit to three of the objects, and a fair fit to the fourth object, 
an ultraluminous infrared galaxy. We find that low incident ionising fluxes ($ H < 10 18 cm~ 2 s _1 ), and 
high gas densities (n H > 10 12 cm~ 3 ) are required to reproduce the observed hydrogen emission line 
ratios. This analysis demonstrates that the use of composite spectra in photoionisation modelling is 
inappropriate; models must be fitted to the individual spectra of quasars. 
Subject headings: galaxies: active - infrared: galaxies - quasars: emission lines 



1. INTRODUCTION 

Physical models of the broad emission line region 
(BELR) of active galactic nuclei (AGNs) have been frus- 
trated by the complexity of the expected emission mech- 
anisms, and a poor understanding of the physical state, 
dynamics and spatial distribution of the line emitting 
gas. Innovativ e new techniqu es, including reverbera- 
tion mapping (TPetersonl Il994f) and microlensing (e.g. 
lO'Dowd et al.ll2011t iBate et all 120081) . are providing in- 
direct methods to measure the physical sizes of the emis- 
sion regions. To date, each of these techniques has only 
been applied to a select sample of sources. However, 
useful models of the BELR will require constraints on all 
physical parameters. Broad emission line fluxes provide 
a powerful diagnostic of the physical conditions in the 
BELR, particularly when coupled with photoionisation 
simulations from codes, such as Cloudy. 

Broad emission lines are a prominent feature in lu- 
minous AGNs, and a defining characteristic of quasars. 
These emission lines are typically broadened by gas ve- 
locities of ~10 4 kms _1 , and are understood to arise from 
gas close to the central black hole and its accretion disk. 
The broad emission line gas reprocesses a substantial 
amount of the light emitted from the central region. The 
emitted line fluxes, and in particular their ratios, are 
strongly dependent on the local physical conditions. The 
broad emission line region is the closest material that 
interacts strongly with the accretion disk that can be 
probed observationally, and a detailed physical under- 
standing is essential in understanding the role of AGN 
energy injection at larger scales. 

A compact broad emission line re gion has been inferred 
from reverberation m apping (e.g. . iBlandford fc McKed 
Il982t iPetersonl fl994f) . As a result of the small spatial 
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scale of the BELR and large distances to the nearest 
quasars, direct spatial imaging of the BELR is not pos- 
sible. In the absence of detailed structural information, 
we can develop photoionization models to predict spec- 
tral line ratios and infer the properties of the emitting 
gas. Such calculations are particularly sensit ive to the 
densi t y and ionisation state of the gas (e.g. iDavidsonl 
Il972t [Osterbrock fc Mathewi 119860 . This provides an 
alternative estimate of the physical scale of the BELR. 

Early BELR photoionisation simulations were lim- 
ited to single cloud models with a single number den- 
sity, column density and ionisation parameter. Sim- 
ulating the entire range of possible BELR parameter 
space was not possible, due to computational limita- 
tions. These simulations were reasonably successful in 
reproduci ng a number of observed emission line flux ra- 
tios (e.g. [DayidIon| [1973; E^anXKrohS [l98l|) . How- 
ever, as more spectroscopic data became available, there 
was conflicting evidence: on one hand, the presence of 
semi-forbidden transitions, such as C hi] A1909 suggeste d 
reasonably low number density gas (|Osterbrocklll970D . 
while the hydrogen emission line ratios and strong Fe II 
emission suggested the pre s ence of high number densit y 
gas (e.g. iRees et all fl989l iCollin-Souffrin et al.1 fl982l) . 
Broad emission line-continuum reverberation measure- 
ments find a wide range in emission line lags relative to 
the continuum variations and line widths, arguing for 
a distributed geometry of physica l conditions (e.g., see 
iKrolik et al.l Il99ll [Peterson 20 y]) . To recon cile these 
conflicting observations, iBaldwin et al.l (|1995| ) proposed 
a locally optimally emitting cloud (LOC) model, where 
gas 'clouds' with a range of physical properties exist at 
each radius in a spatially extended BELR. This model 
takes advantage of large grids of photoionisation models 
that can now be simulated over a h uge range of parame - 
ter space in the density-flux plane ([Korista et al.l fl997). 

The LOC model has been successful in producing ul- 
traviolet and optical emission line flux ratios that match 
observations (IBaldwin et al.ll 19951 iKorista fc Goadll2000l : 
iBottorff et al.l[2000T ) , while maintaining consistency with 
reverberation mapping results (e.g. IKorista fc Goadl 
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2004). The LOC also provides an explanation for the re- 
markably consistent BELR flux ratios, over a very large 
range in AGN luminosity. The name, 'LOC originated 
because many of the optical and ultraviolet (UV) broad 
emission lines emit 'optimally' in one region in parameter 
space. Specifically, the typical emission line is efficiently 
emitted (relative to the incident continuum) over a rela- 
tively narrow region within t he density-incident c ontin- 
uum flux plane (see Fig. 1 of lBaldwin et al.|[l995f ). This 
natural selection effect is due to a combination of ioniza- 
tion potential, collisional de-excitation of the upper level, 
and thermalization at large optical depths. Reverbera- 
tion mapping measurements show that individual emis- 
sion lines respond to variations in the continuum with 
different time delays, and the characteristic radius where 
each line is produced depends on the ionisation potential 
(i.e. high ionisation lines are emitted closer to the con- 
tinuum source). The term, 'LOC,' is less appropriate for 
the Hi and Hei lines, as the production of these lines 
from the continuum rem ains efficient oyer a l arge range 
of the density- flux plane (jKorista et al.l ll997: Ru ff et all 
[20127 in vrevi ). 

The hydrogen and helium lines are of particular inter- 
est because the physical processes governing their pro- 
duction is relatively simple, and there is a negligible de- 
pendence on several free parameters of the model: metal- 
licity, column density and slope of the ionising contin- 
uum. Despite the relative simplicity of the H I emission 
spectrum, some unexplain ed observations persist, such 
as th e Lya/H/3 problem (|Baldwinl 119771 : iNetzer et all 
1995) and the small scatter in measured B almer decre- 
ment s (Ha/H/3) of AGNs with blue continua ([Dong et al.1 
2008). The simple recombination spectrum predictions 
of Case B are invalid for several reasons under BELR 
conditions. Firstly, given the high gas number density, 
high incident flux environment, the optical depths in the 
Balmer and Paschen lines are finite and often large. Sec- 
ondly, the effects of collisions are expected to be very 
important, especially in light of the enhanced lifetime of 
the n = 2 state due to the enormous Lya optical depths 
that are encountered. 

Although the LOC model has been shown to produce 
emission line flux ratios that are consistent with mea- 
sured emis sion lines for optica l and UV metal emis- 
sion lines ([Baldwin et al.l 119951: IKorista fc Goadl 120001) 
as well as helium em ission lines (|Bottorff et al.l l2002t 
IKorista fc Goadl [20041 ) . no comprehensive study of the 
infrared hydrogen lines has been undertaken. The near- 
infrared (NIR) lines, for example the hydrogen Paschen 
series, have the advantage of being affected by dust 
equally, whereas there is differential extinction across 
the optical and UV lines. With advances in infrared 
spectr oscopy, larger sam ples of both low luminosity 
AGN (lLandt et , aLll2008| ) and higher luminosity quasars 
(Glikma ri^t~a^[200fih are being observed. The Paschen 
lines are becoming more widely use d for AGN measu re- 
ments, such as black hole mass (e.g. lKim et"aIll20To[ ). 

This paper investigates a new diagnostic that tests 
the LOC model using near-infrared hydrogen emission 
lines in quasars. The narrow component has a larger 
contribution to the total line flux in lower luminosity 
AGNs, and the systematic uncertainties resulting from 
this subtraction must be quantified. We have concen- 
trated on quasars (as opposed to Seyferts and other 



lower- luminosity AGN), so that deconvolution of the 
emission lines into broad and narrow components is not 
necessary. We use only spectra of individual objects, 
and demonstrate that composite spectra are not suitable 
for determining the physical conditions in quasars. Sec- 
tion [2] describes the photoionisation simulations and the 
simple LOC model used in this analysis. Section [3] de- 
scribes data used. Measurements are compared to the 
LOC model in Section [4j A discussion of the results is 
given in Section [S[ before a brief summary in Section [5] 

2. PHOTOIONIZATION MODELS 



Fol l owing the stan d ard L OC approach ([Baldwin et all 
[19951: IKorista et all [19971 IKorista fc Goadl I2000D . a 
large grid of photoionisation models was computed us- 
ing the photoion isation code, Cloud y version 08.00, 
last described by iFerland et all (|1998l ). Version 08.00 
of Cloudy implemented full resolution of quantum l- 
le vels 5 for the H and He-like sequences for the first 
time (IFerguson fc Ferlandl [19971 IFerguson et al.l 120011 : 
IPorter et al.1 120051 : IPorter fc Ferlandl I2007D . greatly im- 
proving the accuracy of the predicted emission. For 
greater predictive accuracy, the number of resolved quan- 
tum levels in the H I and He I atoms was increased from 
the Cloudy defaults 6 . A brief outline of the grid simu- 
lations is given below. A more detailed des cription of 
the photoionis ation simulations is given by IRuff et al.l 
([2012. in vrev\ ). 

In the context of an LOC model, the term, "cloud" 
refers to an individual emitting body of gas within 
the BELR. Emission line profiles are observed to be 
very smooth, imp lying either a large number of clouds 
(jArav et al.1 119981 ). or a smoothing mechanism. Micro- 
turbulence of ^100 km s -1 would broaden the local line 
profile sufficiently to smooth the bulk lin e profile to 
match the observed ( e.g., as men t ioned in lArav et al" " 



19981: see also iBottorff et all 120001 : iKollatschnv fc Zetz] 



20111) . Micro-t urbulence also repro duces the observed 



UV Fen bump ([Baldwin et al.ll2004D . The presence of a 
small amount of electron scattering (e.g. lEmmering et al.1 
Il992t lLaorll2006l ) or macro-turbulence in the cloud veloc- 
ity field could also smooth the profile. The precise na- 
ture of the line-emitting entities is not the focus of this 
investigation. Instead we focus on constraining general 
physical conditions of the emitting gas using hydrogen 
emission line flux ratios. 

2.1. Grid simulations in the density-flux plane 

Models that predict hydrogen and helium emission line 
strengths were generated for a range of hydrogen num- 
ber density (n H ) and hydrogen ionizing flux ($h) values. 
The plane formed by these parameters represents the 
range of cloud densities and distances from the ionizing 
continuum that are expected to exist within the BELR. 
The output emission of the photoionisation simulations 
is computed as an equivalent width in angstroms relative 
to the incident continuum at 1216A (W\2w)- The equiv- 
alent width is a useful description of the line emission, 
as it measures how efficiently the line is produced from 
the incident continuum. 

5 The azimuthal quantum number, I, is one of the numbers that 
describes the quantum state of an electron. 



6 Simulations used H I and He I model atoms 
resolved levels, respectively. 
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For each point in the grid, Clou dy's standard AGN 
ionisin g continuum (similar to the Ma thews fc Ferlandl 
(1987) continuum) was used. Despite a range in observed 
ionising continua, the H and He emission lines are rela- 
tively insensitive to the shape of the ionising continu um, 
unlike the metal emission lines (jKorista et al.lfl997T ). A 
constant number density within each cloud was assumed 
and a constant column density of iV H = 10 23 cm" 2 was 
also assumed. Although a range of column densities is ex- 
pected within the BELR, the hydrogen and helium emis- 
sion lines are not very sensitive to the clou d column den- 
sity in the range 10 22 < 7V H (cm ~ 2 ) < 10 24 (jKorista et al.1 
Il997t iRuff et al.ll2012. in prepi ). Note that solar chemi- 
cal abundances were used for all photoionisation simula- 
tions. 

2.2. Limits on the number density, ionising flux and 
ionisation parameter 

The range of values for each of the physical param- 
eters in the LOC model is motivated by observations 
and physical considerations. Broad, semi-forbidden C III] 
A 19 09 is observed, implying t he presence of ~10 9 cm~ 3 
gas (|Davidson fc Netzenll979[ ). as it is collisionally sup- 
pressed at higher densities. However, broad forbidden 
lines are not observed, indicating that the electron den- 
sity in the BELR is higher than the critical density of 
the forbidden lines. This gives a lower limit on the 
number density within the BELR and only clouds with 
n H > 10 7 cm" 3 were considered. While none of the broad 
emission lines can be used to determine an upper limit on 
the number density, analyses of the Balmer lines, whose 
ratios are sensitive to high density gas, and strong Fe n 
emission indicate that high density (~10 12 cm" 3 ) gas is 
also p resent within the BELR (e.g.. lCollin-Souffrin et al.l 
Il982t iRees et al.|[l989ft . At these high densities, colli- 
sional excitation becomes important. The collisional ex- 
citation depends on density, temperature and ionisation 
state in a funda mentally different m anner to the recombi- 
nation Case B ([Ferland et al.|[2009t) . Therefore detailed 
photoionization simulations are required. 

Reverberation mapping measurements show that 
broad emission lines origina te at a range o f radii from 
the central ionising source ((Peterson! Il994j ) . At large 
distances from the continuum source, the ionising flux 
is low and the temperature decreases. Below temper- 
atures of ~1500K dust grains condense, ma king the 
line emissivity low (e.g. iNetzer fc Laorl 1 1993T) . This 
temperature corresp onds to a hydrogen i onising flux of 
$ H ~ 10 18 cm" 2 s" 1 (|Nenkova et all 120081 after substi- 
tuting a fiducial bolometric luminosity). The Hi and 
He I emissivities are strongly dependent on the incident 
ionising flux, unl i ke mo st of the prominent metal lines 
(jKorista fc Goadl [2004) . This is because the optical 
depth of the Hi and He I excited state s increases with 
ionising flux (jFerland et al.lll979L fl992T ). until the cloud 
is fully ionised. Close to the ionising source, the hydro- 
gen and helium atoms will be fully ionised and will not 
produce line emission. As the H I and He I line emissiv- 
ity is low at large $ H values, the upper limit on $ H is 
considered to be unimportant. 

Considering the limits discussed above, a large range 
in BELR parameter space was simulated, spanning seven 
decades in both n H and $ H : 10 7 < n H (cm" 3 ) < 10 14 and 
10 17 < $ H (cm" 2 s" 1 ) < 10 24 , each stepped in 0.25 decade 



intervals. Each model was calculated separately by 
Cloudy, giving a total of 841 photoionization calcula- 
tions. Contour plots of each line, showing th e Wme as 
a function of n H and <£> H will be presented in Ru ff et al.l 
(|2012. in prepj ). A lower limit on the ionisation param- 
eter, 



was also considered. This is motivated by observations 
of the calcium infrared triplet, CaiiXYZ AAA8498, 8542, 
8662. If clouds with very low values of the ionisation 
parameter (log U < —6) existed in the BELR , stronger 
CaiiXYZ (jJolvl 119891 : jFerland fc Perssonlll989h and Nai 
A589 5 emission would be observed (see iKorista et all 
1997. for a range of W\2\% distributions in the density- 
flux plane) . In addition, the Fe n emission would be to- 
tally dominated by resonance transiti ons clustered near 
2400 A and 2600 A (jVerner et al.lll999D . which is dissimi- 
lar to the pseudo-continuum observed. A physical inter- 
pretation of the lower limit on the ionisation parameter, 
fmin is that very dense clouds are not found far from 
the ionising continuum. Although this limit may not ex- 
actly represent the true physical cloud distribution, it is 
a convenient representation in the LOC model. 

These considerations provide only weak constraints on 
the upper n H limit, and lower $ H and U m i n limits in the 
BELR. However, they are very important because of the 
large emissivity of H I and He I close to these limits. The 
H I and He I line ratios are particularly sensitive to the 
lower $ H and upper n H values ($ m i n and n m ax, respec- 
tively). Thus rather than using fiducial values, a range 
of upper and lower limits on $ H and n H was considered. 

2.3. Calculating the total equivalent width 

Following the standard LOC prescription, the flux from 
each point in the simulated parameter space was summed 
to calculate t he emission from eac h species of interest. 
As shown by iBaldwin et al.l (j!995l ) , the line luminosity 
is given by: 

/"^max /'"■max 

ilinc oc / / W 12 i(>{r,n H )f(r)g(n H )dn H dr, (2) 

J r m i n J Tlxnin 

where Wi2w(r, n H ) is the equivalent width of a partic- 
ular line (relative to the continuum at 1216A) from a 
single cloud at radius r, with number density n H - The 
cloud covering fractions as functions of radius and num- 
ber density are given by f(r) and g(n n ), respectively. 
For a mo re complete descripti on of the LOC param- 
eters, see Bottorff et al.l ((2002). In the simplest case 
of LOC integration, the coveri ng fractions are given by 
f(r) oc r" 1 and g(n H ) oc n" 1 (jBaldwin et al.l [19951 : see 
also lMatsuoka et al.ll2007t) . The line luminosity can then 
be simplified to: 

/'log'I'max /*logn ma x 

Li ino oc / / Wi2 W ($ H ,n H )d log n H dlog$ H , 

J log ^min J log n m iri 

where the integration limits are free parameters of the 
model and $ H oc r" 2 . This is simply a sum over each 
point in parameter space between the integration lim- 
its, since the simulated photoionized clouds in the log- 
arithmic grid are weighted evenly per decade. We note 
that this weighting scheme maximizes the distribution of 
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Pa 7 /Pa /S 

He I 1.083ft/Pa /S 

H a/U (3 

He I 1.083^/He I 5876 



*m« = 2300 

* mln = 18.00 

n m „ = 13.00 

n mln = 8.00 
U„,„ = -5.5 



-6.5 -6.0 -5.5 -5.0 -4.5 
log 

Fig. 1. — Flux ratios are shown as a function of varying three 
integration limits: $> m i ai n m ax and U m i a . Each panel shows the 
variation of a different limit, while the other limits are held at 
constant at fiducial values, given in the legend. 

cloud characteristics in gas density and incident ioniz- 
ing ffux - a key characteristic of the LOC model. Power 
law indices in the density and distance weighting func- 
tions (f(r),g(n a )) which deviate substantially from -1 
will take on characteristics of a model with a single gas 
density and ionizing continuum source distance 7 . 

This sum (Equation 13]) then allows us to compute the 
predicted total equivalent width (Www) f° r each emis- 
sion line, assuming full geometric coverage. The model 
integrated covering fraction is then determined from the 
ratio between the observed W^i2i6(Lya) and its predicted 
value, above. While in the present work we are mainly 
interested in comparing hydrogen Balmer and Paschen 
emission line flux ratios (which are independent of the in- 
tegrated covering fraction) , Table [5] lists the value of the 
integrated covering fraction that results in a predicted 
equivalent width of Lya of 100A, a typical value (see 
also Section 03]). 



2.4. Ratios calculated from varying the integration 

limits 

Broad emission line flux ratios were predicted for sev- 
eral emission line pairs using Equation |3] Rather than 
using fiducial values of the integration limits to predict 
ratios, a range of upper and lower bounds on each of the 
free parameters: $ m in, fimax and U min was considered. 
From the discussion in Section \2.2\ the $ m i n , "max and 
C^min integration limits were varied over the ranges: 



log$ n 
log n m 
log{7„ 



= [17.00,20.00] 
= [10.25,14.00] 
= [-6.5,-4.51 



(4) 



7 In studies 
IjKorista 



of t he UV emission lin es in NGC 5548 
Goad 2000) and SDSS quasars UNagao et al.l I200r3) 



there is evidence for favoring adoption of a slightly steeper power 
law index in f(r), which tends to favor gas at smaller distances 
from the central continuum source. Here, however, given the lim- 
ited number of observational constraints (discussed in Section 
we choose to adopt the simplest weighting scheme in the LOC in- 
tegration. 



To limit the number of assumptions in the calcula- 
tion, these integration limits span a larger range than 
suggested by the observational limits discussed in Sec- 
tion 12.21 The $min and n max integration limits were 
stepped in 0.25 decade intervals and the U m i n limit in 
0.5 decade intervals. The $ max and n m ; n limits were also 
varied, however the effect on the total emission was small, 
because the H I and He I lines have low emissivity in this 
region of parameter space. The values of $ ma x and n m i n 
were fixed at 10 23 cm _2 s _1 and 10 8 cm -3 , respectively. 
The change in predicted ratio caused by varying these 
limits was less than 1%. 

To visualise how the ratios change as the $ m in, 
"■max and U m - m integration limits are varied, several 
flux ratios are plotted against each of the integra- 
tion limits. Figure [T] plots several predicted ratios: 
Pa a/Pa /3, Pa 7 /Pa /?, Hei 1.083/x/Pa /3, Ha/H/3 and He i 



1.083/x/Hei A5876 as a function of $ n 



and U n 



integration limits. The three panels show ratios plotted 
as a function of each of the integration limits over the 
range given in Equation 0] while the other integration 
limits are held constant at fiducial values (given in the 
legend). 

3. APPLICATION TO DATA 

Very few quasars have b een observed in the in frared 
with high signal-to-noise. iGlikman et al.l ([2006) pub- 
lished a set of 27 low signal-to-noise quasar spectra, that 
were presented as a composite quasar spectrum. Aver- 
aging quasar spectra is complicated by several factors, 
which limits the accuracy of emission line fluxes mea- 
sured from the composite spectrum, see Section 15.31 for 
details. 

Rather than using the composite spectrum, individ- 
ual s pectra with the h ighest signal-to-noise ratios (SNR) 
from IGlikman et al.l (|2006| ) were analysed. In prac- 
tise our cutoff was SNR > 10 per pixel in the contin- 
uum between the Pa a and Pa/3 emission lines (exclud- 
ing any regions of significant atmospheric absorption). 
Although there have been other NIR observations of 
Seyferts (jLandt et al.l 12008') . only quasars were studied, 
as the subtraction of a significant narrow component in 
Seyfert galaxies adds additional uncertainty, which must 
be quantified in order to compare measured to predicted 
line fluxes. 

3.1. Measured broad emission line fluxes 

Four quasar spectra from the IGlikman et al.l {2006) 
sample have high enough signal-to-noise ratios to mea- 
sure the strongest Paschen emis sion line fluxes. For each 
emission line in the SDSS and IGlikman et al.l spectra, 
the local continuum was measured over a wavelength 
range ±200 — 500A, centred on the line's expected wave- 
length. Linear fits to the continuum were performed after 
removing the line itself and any 2a outliers. To measure 
line strengths, both gaussian fits and simple integrations 
of the flux per wavelength bin were performed for each 
line. While gaussians are an imperfect match to the line 
shapes, they minimise parameterisation while allowing 
for de-blending of the Hei and Pa 7 lines. Integrating 
the flux per wavelength bin over ±2 FWHM reproduced 
the gaussian line fluxes within the errors. Summing over 
±3 FWHM admits smaller outlying features and so tends 
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Fig. 2. — Glikman NIR spectra (dark grey), scaled SDSS spectra (black), GALEX far and near UV measurements (error bars). The 
extrapolated optical continuum (dotted line) is also shown, with an estimate of the 1216A continuum flux marked with a diamond. The 
insets show closeups of the Pa7+Hcl 1.083/i, Pa/3 and Paa emission lines with best fit gaussians. 



TABLE 1 

Measured NIR and optical line fluxes 



Object 



He i 5876A 
F e(F) 



H7 434lA 
F e(F) 



H/3 4861A 
F e(F) 



Ha 6563A 
F e(F) 



He 1 1.083/xm 

F e(F) 



Pa 7 1.094/Ltm 
F e(F) 



Pa p 1.282/im 
F e(F) 



Paa 1.876/im 
F e(F) 



J005812+160201 
J010226-003904 
J015950+002340 a 
J032213+005513 



347 10% 1285 10% 4364 10% 11446 10% 1580 14% 626 17% 1078 14% 1388 12% 

630 11% 2577 10% 6281 10% 18003 10% 2776 11% 1264 10% 1370 10% 1661 10% 

2452 10% 4259 10% 7506 10% 13447 50% 5222 10% 1296 11% 2100 11% 2772 11% 

682 10% 3883 10% 7888 10% 22781 10% 2202 11% 1151 11% 1779 10% 2421 11% 



NOTE. Measured NIR and optical line fluxes (10~ 17 ergs -1 cm" 2 ) with their associated measurement uncertainties. The optical lines: Ha, H/3, H7 and He 1 A5876 were 

measured directly from the SDSS spectra, and have been scaled by fgDSS from Tablc[2]for comparison with the infrared lines. A 10% systematic uncertainty was estimated on each 
emission line flux, predominantly due to use of Gaussian fitting (see text). JO 15950 + 002340 Q is a ULIRG with an extremely asymmetric Ha lino; the 50% error on Ha accounts for 
the inadequacy of a Gaussian fit for this line. 



to result in an over-estimate of the line strength. Mea- 
sured line fluxes are shown in Tabic [1] together with their 
1-ct measurement uncertainties. The SDSS line strengths 
were scaled by /sdss to allow direct comparison. An ad- 
ditional 10% systematic error was conservatively added 
in quadrature to measurement errors in all figures and 
calculations. 

3.2. Measured continuum fluxes 

For each of the four objects in the sample, SDSS spec- 
tra and GALEX ultraviolet flux measurements were ob- 
tain ed from archival dat a. Colour excess values, E(B — 
V) (|Schlegel et alJll998t l are given in Table and these 
were use d to correct all the data for Galactic extinction 
following lFitzpatricld (|1999f l . Note that although the con- 



tinuum fluxes inferred in this section have large uncer- 
tainties, the inferred physical conditions do not depend 
on these values. The measured continuum flux values 
are only used in Section 14.41 to check that the model 
produces enough energy to be consistent with previous 
observations. 

For each object, the continuum level of the SDSS spec- 
trum was scaled to match the NIR continuum across the 
region of overlap. The scaling factor for each SDSS spec- 
trum, denoted /sdss, was typically measured to better 
than 2% accuracy. The continuum flux at 5100A ($5100) 
was measured from the SDSS spectrum, which can be 
scaled to match the NIR data using /sdss- An extrap- 
olation of the best linear fit to the optical continuum 
was made to give an estimate of the rest-frame 1216A 
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TABLE 2 

OBSEHVABLES, FLUX-MATCHING AND ESTIMATES OF THE 1216A AND 5100A CONTINUUM FLUXES 









GALEX 








SDSS 




Object 


z 


Sfuv 


Snuv 


Eb-v 


5'l216 


S1216 


S5100 


/SDSS 


^(/sdss) 


SDSSJ005812. 8+160201. 3 


0.211 




20.6 


0.075 




110 


42 


0.729 


1.2% 


SDSSJ010226. 3-003904. 6 


0.295 


880 


431 


0.036 


2300 


2200 


130 


1.452 


0.6% 


SDSSJ015950. 2+002340. 8 


0.163 


1020 


580 


0.029 


2500 


2600 


170 


1.073 


0.6% 


SDSSJ032213. 8+005513. 4 


0.185 


164 


108 


0.120 


790 


1100 


100 


1.003 


0.8% 



Note. — SDSS target name; redshift; GALEX raw observed fluxes (10 17 ergs 1 cm 2 A x ) at 
1516 ± 134A (FUV) and 2267 ± 366A (NUV); Galactic extinction color excess ijSchlegel et al.|[l9M ): 
GALEX estimate of continuum flux at 1216A; SDSS-extrapolated continuum flux at 1216A and measured 
continuum at 5100A; flux factor to match SDSS to the NIR spectra and its uncertainty, e. Note that 
GALEX Sprjv an( i ■S'nuv fluxes have not been corrected for extinction in this table. 



,10058 12+ 160201 JO 10226-003904 
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Fig. 3. — Contours of reduced \ 2 in the <I> m i n — Hmax plane for each 
of the four objects. Measured Pa a/Pa /3, Pa7/Pa/3 and Ha/Pa/3 
ratios were used to calculate \ 2 ■ I n each panel, the location of 
the minimum \ 2 la shown with a red diamond (see Table \3\ for 
minimum values of reduced X 2 )- 

continuum flux (5i2i6(SDSS)). <Si2i6(SDSS) is likeiy to 
overestimate the true continuum flux at 12 16 A due to the 
turnover in the spectrum. An alternative estimate of the 
1216A continuum flux was obtained from the GALEX 
far-ultraviolet (FUV) measurement. Due to the shape, 
breadth, and variable position of Lya in the GALEX 
bandpass for our sample, we cannot use the equivalent 
width of Lya alone to estimate the 12 16 A continuum 
flux. Instead we convolve a n appropriately redshifted 
composite quasar spectrum ()Vanden Berk et "aLl 120011) 
with the GALEX FUV bandpass to estimate the con- 
tinuum level. The SDSS and GALEX estimates of the 
1216A continuum flux level are in good agreement, con- 
sidering the significant uncertainty in extrapolating the 
continuum. Figure [2] shows the NIR spectra, scaled SDSS 
spectra and GALEX photometry. 

4. RESULTS 

The primary aim of this paper is to find a set of integra- 
tion limits for the hydrogen number density and hydro- 
gen ionising flux (see Equation [3]) that best fit the mea- 



sured infrared hydrogen emission line ratios for individ- 
ual quasars. Since there are only a few measured infrared 
ratios, only a small number of model free parameters can 
be fitted. Thus the aim of this analysis is to illustrate 
the potential of using optical and infrared hydrogen and 
helium emission lines as probes of the physical conditions 
within the BELR, rather than constraining a definitive 
model. 

4.1. Best-fit integration limits 

A set of 'best-fit' integration limits was found for each 
object using the measured Pa a/Pa/?, Pa7/Pa/3 and 
Ha/Pa fj ratios. Ratios involving optical emission lines 
were not used, as they suffer differential extinction from 
dust and the amount of dust is unknown. There are only 
three independent infrared hydrogen ratios, so the model 
must be limited to fitting to two parameters. Since the 
ratios are fairly constant with the lower limit on the ioni- 
sation parameter, log t/ m m = —5.5 is fixed, and and 
n max are free parameters. 

The measured emission line ratios of four objects: 
J005812+160201, J010226-003904, J015950+002340 
and J032213+005513 were used to find the most likely 
set of $ m in and n max limits by minimising x 2 : 

2 \ ^ (-^obs,i -^model,z) /_\ 

* = 2^ — —2 , (5) 

stat,z 1 sys,z 

where i? bs arid i? mo doi are the observed and predicted 
ratios, and t7 s tat and cr S y S are the statistical and system- 
atic uncertainties on each measured ratio. The sum is 
over each ratio used to find the best model: Pa a/Pa f3, 
Pa 7 /Pa j3 and Ha/Pa (3. R mo dc\ is the predicted ratio 
for each combination of $ mm and n max values calculated 
using Equation [31 over the ranges given in Equation 2J 
Table [3J shows the minimum \ 2 value for each object. 
Figure [3] shows the \ 2 contours for each object in the 
^mm-imax plane. This plot shows that the % 2 value is 
minimised where the density upper integration limit is 
high and the lower limit on the ionising flux is low. 

Figure |4] shows the best-fit simulated ratios as a func- 
tion of the two free parameters: rt max and $ m i n . The 
best-fit n max and 3> m in limits are given in the legend for 
each object. The n m ; n and f max integration limits were 
fixed at fiducial values given in Section 12.41 The mea- 
sured infrared ratios are over-plotted as a band, indicat- 
ing a lfi uncertainty. Note that the Hei 1.083/j,/Pa/3 
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ratio was not used as a constraint, as this ratio predicts 
$ m in and n max values that are systematically offset from 
those given by the hydrogen line ratios, as shown in Fig- 
ure m 

4.2. Predicted Balmer ratios 

The limits on n max and $ m i n that were calculated by 
analysing the NIR emission lines were then used to pre- 
dict Ha/H/3 and H7/H/3 ratios. The predicted ratios 
are shown with the measur ed ratios in Table El Two 
of the four objects are in the lDong et al.l f|2008T ) sample, 
that selected AGNs with blue optical continua in order 
to minimise the effect of dust extinction. If blue con- 
tinua do indicate less dust along the line of sight, the 
simulated ratios should be in good agreement with the 
measured ratios for these objects. This is true for the 
two objects listed, however, a much larger sample with 
a larger number of ratios studied is required to test this 
hypothesis. J015950+002340 also has a blue continuum, 
however simulated values cannot be compared to mea- 
surements, due to the poor fit to the Ha line. 

For all objects, the predicted Balmer decrements are 
slightly larger than the measured values, which is the op- 
posite effect to dust. However, the predicted values are 
within measurement uncertainties, which may indicate 
that there is little internal reddening in these objects. 
The predicted H7/H/3 ratios are consistent with the mea- 
sured data. The results for each object are summarised 
briefly in Section l4~5l 

4.3. The predicted Lya/H/3 and UV line ratios 

Although UV emission lines have not been measured 
for these objects, a selection of predicted UV emission 
line ratios are presented in this section. Typical mea- 
sured strengths of stronger UV lines were compared to 
the predicted values to check the consistency of our spec- 
tral model derived from the measured hydrogen spec- 
trum. We note that the UV metal emission line strengths 
depend more on the shape of the incident ionising con- 
tinuum (and to varying degrees the gas metallicity) than 
do the B almer and Paschen hydrogen lines under study 
(e.g., see lKorista et al.lll998l for a discussion). The pre- 
dicted UV emission line strengths reported here are those 
that result from the best-fit in the physical parameters 
(n max , $ m in) determined from the three hydrogen emis- 
sion line ratios (as noted in Table [3]) in each of the four 
AGN. 

Predicted Lya/H/? ratios for each object are listed 
in Table |4j Most of the Lya/H/3 predicted ratios are 
lowe r than the simple, classical photoionisation predic- 
tion (jOsterbrock fe Ferl and 2006). This is because of the 
inclusion of high density gas. Lya thermalises at high 
densities, whereas the Balmer lines continue to emit effi- 
ciently, causing the predicted Lya/H/? ratio to decrease 
with the inclusion of high n m gas. Nevertheless, these ra- 
tios remain larger than many of those measured in AGN, 
and this remains an unsolved problem. However, uncor- 
rected extinction due to grains extrinsic to the Milky 
Wa y may explain muc h of the remaining discrepancies 
(see iNetzer et aT1ll995f ) . 

A selection of predicted UV metal lines (relative to 
Lya) are listed in Table SJ These lines were chosen as 
they are relatively strong, and they span a large range 



in ionisation parameter. The Ovi, Civ and Mgn emis- 
sion lines are also the most important coolants within 
their ionisation zones (excepting emission by Fen). As 
these lines are major coolants, the line strengths are 
not very sensitive to the gas metallicity, due to the 
thermostatic effects of such lines. The predicted C iv 
values a re somewhat highe r than the ratio measured 
from the lVanden Berk et al.l (|2001h composite spectrum, 
however, the predicted Civ/Lya ratio is strongly de- 
pendent on the hardness of the ionising continuum. 
The observed value shows significant varia tion; for ex- 
ample , the Civ/Lya ratio is 0.62 in the IZheng et al.1 
(11997ft composite, and is c loser to 1 in Seyfert 1 spec - 
tra (iKorista fc Goadl [2000T; see also lOsmer et al.lll994l) . 
Although additional free parameters, such as metallicity 
and the shape of the ionising continuum, have not been 
considered here, the predicted strengths of the UV metal 
lines relative to Lya are generally consistent with typical 
measured values. 

4.4. Do the models produce enough energy? 

A simple check, albeit with a substantial uncertainty, 
can establish whether the best-fit models produce enough 
energy to reproduce the observed total integrated line 
flux. The largest uncertainty in comparing measured 
equivalent width values to predicted line strengths is 
knowing which continuum flux to compare the line flux 
to. The NIR continuum is a poor indicator of the inci- 
dent ionising continuum (i.e. photons with wavelengths 
shorter than 912A), and is contaminated with hot ther- 
mal dust emission as well as that of the host galaxy. To 
test whether the models produce enough energy, the pre- 
dicted and measured Pa (3 equivalent widths were com- 
pared, relative to two different points in the continuum: 
1216A and 5100A. Note that the predicted equivalent 
width depends on the integrated covering factor. In this 
section, the integrated covering factor was calculated, 
such that the simulated equivalent width of Lya is 100A. 
This is a typical measured value for IVi2i6(Lya), and was 
chosen to ensure that the model produces enough energy 
to be consistent with the observed value. For each ob- 
ject, this value of f c was used to calculate all equivalent 
width values presented in this section. 

The equivalent width relative to the continuum flux at 
1216A is often chosen because it is easily observable and 
close to the ionising continuum. Contaminants such as 
light from the host galaxy and dusty torus are negligible 
in this region. The measured value was estimated using 
both the GALEX and SDSS UV continuum flux values 
(51216) listed in Table [2] The observed rest- frame equiv- 
alent width of Pa f3 relative to the incident continuum at 
1216A (in A) is given by 

' (6) 



Wm o*n = x 



#1216 (1 + z) 

where 5*1216 is the measured incident continuum flux at 
1216A and F(P&/3) is the total integrated emission line 
flux, given in Table [T] The dominant measurement un- 
certainties are quasar variability and scaling the UV mea- 
surement to the NIR. There is an additional uncertainty 
from dust attenuation. The measured 5*1216 values show 
a large variance between the GALEX and SDSS esti- 
mates, therefore this value is estimated to be accurate 
within a factor of 2. 
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Fig. 4. — Measured ratios are plotted compared to the best-fit models. The solid, dashed, dot dashed and dotted lines show the simulated 
Paa/Pa/3, Pa7/Pa/3, Hel 1.083/i/Pa/3 and Ha/Pa/3 ratios, respectively. The measured values are shown as a band with lc uncertainties 
in a corresponding colour. For legibility, the very large Ha measurement uncertainty is not shown for J015950+002340, although it was 
used in the best-fit calculation. The parameters that produce the best fit for each object have been plotted and the values are given in the 
legend for each object. Note that a log scale has been used on the y-axis, so that the range of ratios can easily be seen. 



The measured Pa/3 equivafent width vaiue reiative to 
the incident continuum at 5100A (VT 5 ioo(Pa/3)) was aiso 
calculated. At this wavelength the continuum flux is di- 
rectly measurable from the scaled SDSS spectrum (unlike 
the SDSS continuum at 1216A 2 which is an extrapolated 
estimate). However, the 5100A continuum is much fur- 



ther from the ionising continuum and host galaxy con- 
tamination is not negligible at this wavelength, making 
it a poorer proxy for the ionising continuum. There may 
also be some thermal emission from the BELR clouds and 
reflected light from the dusty torus contributing to the 
measured continuum flux at this wavelength. The calcu- 
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TABLE 3 
Measured and predicted ratios 











Used for fit 








Not used for fit 




Object 


x 2 


Pa a/Pa 


P 


Pa7/Pa 


P 


Ha/Pa 


P 


Ha/H/3 


H7/H/3 


He 1 10830/Hei 5876 






m 


P 


m 


P 


m 


P 


m p 


m p 


m 


P 


J005812+160201 


1.5 


1.29L0.2 


1.34 


0.58L0.1 


0.74 


10.6L2.2 


10.6 


2.62L0.4 3.04 


0.29±0.1 0.53 


4.54L0.8 


3.09 


J010226-003904 d 


0.27 


1.21±0.2 


1.18 


0.92±0.1 


0.86 


13.1±2.0 


13.4 


2.87L0.4 3.18 


0.41±0.1 0.46 


4.41L0.7 


2.42 


J015950+002340 


3.1 


1.32±0.2 


1.30 


0.62±0.1 


0.74 


6.40±3.0 


10.3 


1.79L0.9 2.90 


0.57±0.1 0.55 


2.13±0.3 


2.96 


J032213+005513 d 


0.90 


1.36±0.2 


1.44 


0.65L0.1 


0.73 


12.8L2.2 


12.8 


2.89±0.4 3.60 


0.49±0.1 0.45 


3.23L0.5 


3.23 



Note. — Measured (m) and predicted (p) broad emission line flux ratios. Using the best-fit parameters discussed in Section 14.11 the 
Ha/H/3, H7/H/3 and He I 1.083/i/Hel A5876 ratios were predicted for each object. The pre dicted ratio s are c ompared to the measured 
SDSS ratios. A 10% systematic was included on each measured flux. Objects that are in the lDong et al.l (|2008T ) sample are denoted by d. 



TABLE 4 

Predicted LYa/H/3 and ultraviolet line ratios 



Object 


Lya/H/3 


O vi+Ly/3/Lya 


Civ/Lya 


Mgn/Lya 


J005812+160201 


19 


0.14 


0.33 


0.22 


J010226-003904 


40 


0.24 


0.40 


0.09 


J015950+002340 


17 


0.14 


0.33 


0.22 


J032213+005513 


31 


0.16 


0.35 


0.16 


measured" 


12 


0.10 


0.25 


0.14 



Note. — Predicted Lya/H/3, and Ovi A1035+Ly/3, Civ A1549 and 
Mgn A2800 line ratios relative to Lya. The line ratios were predicted 

using the best-fit parameters derived i n Section 14.11 

a The measured line ratios are from the lVanden Berk et all i|200lT ) com- 
posite spectrum. 



lation of Wsioo(Pa/3) is analogous to the calculation of 
Wi2i6(Pa/3) given in Equation[B] Given the large uncer- 
tainties in scaling the continuum measurements to the 
NIR, quasar variability and host galaxy contamination 
in the 6*5100 measurement, the measured values listed in 
Table [5] should only be considered accurate to within a 
factor of 2. 

Using the best-fit n ma x and $ m i n integration limits de- 
termined in Section |4~T1 the simulated Wm6 (Pa/?) and 
W5ioo(Pa P) values were calculated for each object. Both 
measures of equivalent width can be output by Cloudy 
directly. (Note that these estimates assume the standard 
Cloudy continuum described in Section \2. 11 ) To compare 
the simulated equivalent width values to the measured 
values, predicted Wm6(Pa {$) and W / 5ioo(Pa/3) values 
must be scaled by the integrated cloud covering fraction 
(/ c ). For each object, the integrated covering fraction 
required to produce VLi2i6(Lya)=100A (a typical value) 
was computed for the best-fit parameters. This value of 
f c was then used to scale the predicted equivalent width 
values. 

The measured and predicted Wi2ie(Pa/3) and 
W / 5ioo(Pa/3) values, and the value of f c are listed in Ta- 
ble [5J As shown in this table, the computed integrated 
covering fraction required to produce Wi2i6(Lya)=100A 
for each object is not unrealistically large. The measured 
and predicted W\2i§ (Pa f3) values are also in good agree- 



ment (with the exception of J005812+160201, which has 
a shallower UV continuum slope than the Cloudy AGN 
continuum). Therefore Table [5] shows that the best- fit 
models produce enough energy to be consistent with both 
the measured Pa ft and previously observed Lya equiva- 
lent width values. 

4.5. Results for each object 

A brief description of the results for each object is given 
below. 

J005812+160201: the model is a good fit to the data, 
although the predicted H7 /H/3 ratio is 2a larger than the 
measured ratio. The measured H7/H/3 ratio is low com- 
pared to the other objects. The predicted total equiva- 
lent width barely produces enough energy to be consis- 
tent with the observations. The low predicted equivalent 
width of Pa j3 is at least partially due the assumed Cloudy 
continuum, which is a poor approximation in this case, 
as the measured UV spectrum is much shallower than 
the Cloudy continuum. The measured W / i2i6(Pa/3) is 
high compared to the other objects, again because the 
continuum slope is very shallow. 

J010226-003904: the model is a good fit to the data. 
This object has a strange BAL-like feature in He 1 1.083/1, 
probably due to the atmosphere. Howeve r, this feature 
could be real, similar to that reported by iLeighlv et al.l 
(|2011l) . The predicted ^min value is quite high, which 
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TABLE 5 

Predicted rest frame equivalent width of Pa/3 



Object 


Wi2i 6 (Pa 


P) 


VK5100 


(Pa/3) 


fc 




GALEX 


SDSS 


P 


SDSS 


P 


p 


J005812+160201 




10 


1.5 


30 


12 


0.31 


J010226-003904 d 


0.3 


0.3 


0.57 


6 


4.5 


0.41 1 


J015950+002340 


0.7 


0.6 


1.7 


10 


13.3 


0.32 


J032213+005513 d 


2 


1 


0.89 


10 


7.1 


0.30 r 



is similar for all objects. High $ m i n limits are gener- 
ally highly disfavoured. This suggests that log$ m ; n < 
18, and that low incident ionising fluxes are required 
: to reproduce observed hydrogen emission line flux ra- 
tios. The determined limits on ^ min are also consisten t 
with the dust sublimation radius (Ncnk ova et alJ feOOS). 
This is in agreement with several recent studies that 
suggest the Balmer lines and Mgli in large part come 
from the ou ter BELR and near or at the dust sublima- 
tion radius (IZhu et al.ll200l ICzernv fc Hrvniewic3l2011l: 



Note. — Measured and predicted (p) Pa/3 equivalent width 
values (in A) normalised to two different points in the continuum: 
1216A (SDSS and GALEX) and 5100A (SDSS only). The pre- 
dicted equivalent width values are calculated using the best-fit 
parameters and were scaled using an integrated covering fraction 
that gives a predicted Lya equivalent width of 100A. The value 
of f c required to produce Wi2i6(Lya) = 100A is also given. Uncer- 
tainties on the measured equivalent width values are large, and 
these values are estimated to be accurate to a factor of ~2. 



results from the large measured Pa 7 /Pa f3 ratio. 

J015950+002340: the model is a fair fit to the data, 
as the reduced % 2 value is 3.1. This source, also known as 
Mrk 1014, is an ultra-lu minous infrared galaxy (ULIRG) 
(see lArmus et al.ll2004l and references therein). The Ha 
emission line is highly asymmetric and is not well fit by 
a Gaussian, making comparisons to simulations difficult. 
The high density predicted could be a consequence of 
the underestimation of the Ha flux. Note that the mea- 
sured Ha/Pa (3 and Ha/H/3 ratios are low compared to 
the other objects. Also note that the optical and NIR 
He I line flux (relative to the hydrogen line flux) is high 
compared to the other objects. 

J032213+005513: the model produced a good fit to 
the emission line ratios, although the predicted Ha/H/3 
ratio is 2a above the measured value. 

5. DISCUSSION 

The predictions of Cloudy LOC models have been 
tested by comparing simulated ratios to measured near- 
infrared hydrogen emission line ratios. Near-infrared hy- 
drogen emission lines were chosen to better understand 
hydrogen emission from the BELR, and because near- 
infrared ratios are insensitive to uncertainties in the line- 
of-sight extinction due to the presence of dust. The 
signal-to-noise of the published spectra is limited, so only 
four quasars were considered suitable. Two questions 
were considered: 

(i) Is a good fit to the observed line ratios possible 
using the LOC model? and 

(ii) What are the integration limits for the best-fit LOC 
models? 

Due to the limited number of independent line ratios, 
only two integration limits in the LOC model could be 
tested. 

5.1. Can a simple LOC model fit these observations? 

In three of the four cases, a good fit to the hydro- 
gen emission lines is obtained by the LOC model. In 
the fourth case the fit is fair. In addition, the shape 
of the x 2 contour distributions in the density-flux plane 



- Figure [3] shows that lower n max values in combination 
with low <E> m i n values are also highly disfavoured (i.e. % 2 
values are high in the lower left corner of the plot), mean- 
ing that BELRs with a low high-density cutoff and low 
low-incident ionising flux cutoff are unlikely solutions. 
The minimum of the y 2 contours tend to lie on a diago- 
nal in the n. max -$ m ; n plane: if n max is larger, then $ m ; n is 
smaller. If the BELR includes regions where the incident 
ionising flux is very low, then gas with a very high num- 
ber density is required to reproduce the observed emis- 
sion line ratios. This solution is quite unlikely though, 
as dust grains condense at low <j> H , destroying the line 
emission. The validity of the inferred limits is uncertain. 
A larger sample, with higher signal-to-noise is required 
to test these results. 

The predicted Ha/H/3 and H7/H/3 ratios are gen- 
erally in very good agreement with the measured ra- 
tios, as shown in Table El However, the measured He 1 
1.083/i/Pa/3 ratio does not match the simple LOC model 
described here. We also find that the predicted Hci 
5876/H/3 line ratio is significantly higher than the mea- 
sured values. The simulated values are typically a fac- 
tor of ~2 greater than the measured values. This offset 
could be caused by the simplicity of the covering frac- 
tions assumed in this model. The shape of the IUi2i6(He 1 
1.083/i) contours is quite different to the Balmer and 
Paschen lines, making ratios involving Hei 1.083/i more 
sensitive to the power laws assumed in Equation [2] Al- 
though the slope of the power laws are not expected to 
deviate significantly from —1, refinement of the models 
using different power law distributions in f(r) and g(n m ) 
will be the subject of future studies. The predicted over- 
estimate of both the He I A5876/H/3 and 1.083/i/Pa/3 line 
ratios could also indicate that Cloudy overpredicts the 
He 1 emission spectrum relative to H. 

5.2. Dust 

The amount of attenuation from internal dust (i.e. 
dust within the host galaxy) is unknown. The Balmer 
decrement is the most widely measured ratio is AGN 
spectroscopy (e.g. IOsterbrocklfl977HGroves et al.ll2012j h 
and has long been used as an in dic ator of the a mount 
of extinction (e.g. iBermanl fl936l ). iDong et all (|2008fl 
measured the Balmer decrement for a sample of 446 
Seyfert Is and QSOs with blue continua, which were se- 
lected to minimise the effect of dust extinction. They 
found only a small intrinsic scatter i n the measu r ed ra - 
tio. The small intrinsic scatter in the IDong et al.l (2008) 
measurement suggests that the intrinsic Balmer decre- 
ment could be relatively consistent between objects, and 
the observed variance is due to different amounts of dust 
along the line-of-sight within the host galaxy. However, 
this result relies on the slope of the continuum being a 
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dust indicator (see lDavis et al.ll2007l for a discussion). 

The continua of our subset of quasars are reasonably 
blue, with the exception of J005812+160201. If dust at- 
tenuation caused reddening of the continuum in this case, 
a larger observed Ha/H/3 would be expected. Unfortu- 
nately it is not possible to test this assertion with such 
a limited sample size. The presence of uncorrected ex- 
tinction due to dust along the line of sight would also 
increase the observed He I 1.083/He I 5876 ratio above 
that of the predicted one (see Table G2) ■ 

With better quality spectra, an estimate of the amount 
of dust attenuation could be made. The Lya/H/3 prob- 
lem could be at least partially expl ained by internal 
extinction (|Shuder fc MacAlpindll979t ) and constraining 
the amount of dust attenuation is important to resolve 
this matter. 

5.3. Composite quasar spectra 

The v alue of the Pa a/P a line ratio measured 
from the iGlikman et al.l (|2006|) composite spectrum was 
0.64±0.01. However, this value is much lower than 
the values measured from individual objects, which are 
consistently above 1.0 (see Table [3]). Other observed 
Pa a/Pa j3 emission l ine ratios in quas ars are consistent 
with this result (e.g. ISoifer et al.l 1200 41 who measured a 
Pa a/Pa /3 ratio of 1.05±0.2 in a z = 3.91 quasar). Ratios 
of Pa a/Pa /3 that are abo ve 1.0 have also bee n measured 
in lower luminosity AGN (jLandt et al.ll2008l ). 

While it is possible that systematic effects cause ob- 
jects with larger Pa a/Pa (3 ratios to be preferentially 
selected in our sample, a bias caused by the compos- 
ite methodology is a more likely explanation. Compos- 
ite spectra are useful for describing generic properties 
of a sample and for identifying emission lines that are 
too weak to be observed in individual spectra. However, 
care must be taken in ascribing physical significance to 
individual features of composite spectra. It is not clear 
whether composite spectra are useful for investigating de- 
tailed physical properties of the line emitting gas. Emis- 
sion line ratios in the composite spectrum can also be 
affected by: poor signal-to-noise ratios of individual spec- 
tra, and wavelength ranges where there a re few objects 
contri buting to the final spectrum. See ISulentic et al.1 
(2002) for a discussion of how emission line ratios are 
affected, and an alternative method for producing com- 
posite spectra. 

5.4. Caveats 

The NIR spectroscopy was not simultaneous with the 
optical and UV measurements. As quasars are intrin- 
sically variable, this will introduce a systematic uncer- 
tainty. Even after continuum matching the line fluxes 
will alter, t his effect is known as th e intrinsic Baldwin 
effect (e.g. iPogge fc Peterson! I1992D . There is often a 
considerable difference between GALEX and SDSS S1216 
values. The objects in Table[3]SDSS are all observed dur- 
ing 2000. NIR observations are from September 2003 to 
September 2004 and GALEX observations are in 2003 
(and one in 2005). 

Fitting the emission lines with Gaussians introduces 
another systematic uncertainty. IGlikman et al.1 ()2006f ) 
point out that Gaussians do not fit the broad wings well, 
however Gaussians generally fit the bulk of the broad 
emission line flux, while requiring only a small number 



of free parameters. The estimated co ntribution fr o m the 
narrow component is less than 5%. iDong et all (2008) 
deconvolved the emission lines into both a broad and nar- 
row component, and in each case the equivalent width of 
the narrow component is only a few percent of the broad 
component. Fe 11 contamination was not considered. Fit- 
ting the Fen contamination would introduce additional 
errors, given the low signal-to-noise of the data. 

Fixing the lower limit on the ionisation parameter, 
f/tnin is not ideal, however, the ratios considered are in- 
sensitive to C/ m i n and with only a few measured ratios, 
the model must be limited in complexity. The l/ m ; n limit 
would be better determined with the Ca 11 infr ared triplet 
(CallXYZ) or other [/ m i„-sensitive ratios (see I Ruff et aTl 
120127 in vrevi for a discussion). Single point solutions 
in the density-flux plane were not considered for reasons 
discussed in Section |2~2"1 

Since the Balmer and Paschen lines are optically thick, 
the geometry will affect the predicted emission. If the 
BELR geometry is such that emission from the 'back' of 
the clouds is not observed, the predicted emission line ra- 
tios will change (iFerland et alj |l992t lO'Brien et"ai] |1994t 
iRuff et alJl20127 in prepi ). 

Any broad emission line flux measurement is from a 
snap-shot in time. The ratios will vary with time as the 
continuum luminosity varies, and this will be dependent 
on the distribution of the gas - both in geometry and 
physical state. Comparing measured flux ratios can be 
problematic if the source is highly variable and lines are 
formed at very different radii. The recombination lines 
are emitted relatively efficiently over a broad range of 
physical conditions, and their continuum reprocessing ef- 
ficiencies are largely dependent on the incident flux, and 
therefore the distance from the ion ising source, due to 
optical depths in the exci ted states ([Korista et alJU997t 
I Ruff et ahl 12012. in prepi ) . This manifests in different 
measured rever beration mapping signals and line profile 
varia bility (e.g. ISergeev et al1l2001t iKollatschnv fc Zetzll 
2010). For example, not only do the Balmer lines show 
different measured time lags, but the line responsivity is 
observed to decrease from line centre to line wing (see 
iKorista fc Goadl [2004 for a discussion). This could be 
corrected for if the distribution of broad emission line 
gas could be determined using high quality spectral and 
tem poral data in conj unction with photoionisation mod- 
els (Home et al.ll2003l ). 

6. SUMMARY 

We have demonstrated a new diagnostic for the den- 
sity and ionisation state of the hydrogen broad line gas 
in quasars. The results should be seen as preliminary, 
being based on fairly low SNR (SNR ~ 10 per pixel in 
the continuum) data for a small number of objects. The 
technique is promising because the hydrogen line ratios 
depend strongly on the physical conditions, but weakly 
on metallicity, column density and slope of the ionis- 
ing continuum. The ionising flux is a good measure of 
distance from the source, and the technique offers an al- 
ternative route to measuring the structure of the quasar 
central engine, in addition to the physical state of the 
broad line gas. 

Despite the simplicity of the model, broad NIR hydro- 
gen emission line ratios can be reproduced by an LOC 
integration. All hydrogen ratios predicted, including the 
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Balmer decrement, were consistent with those observed. 
Low incident ionising flux values ($ H < 10 18 cm~ 2 s ) 
arc required to reproduce the observed hydrogen emis- 
sion line ratios. This predicted lower limit on $ H gives 
an outer radius of the BELR that is cons istent with the 
radius of dust sublimation calculated by Nenkova et al. 
<p008h - 

This model predicts that high number density gas 
(n H > 10 12 cm~ 3 ) is required to reproduce the observed 
hydrogen emission line ratios. The upper limit on the 
gas number density is high compared to previous calcula- 
tions, however, all predictions stated here should be con- 
sidered estimates, because of the simplicity of the model 
and limitations in the quality and number of spectra. 

With a larger sample of high signal-to-noise quasar 
spectra in the NIR, higher order lines in the Paschen 
series could also be observed. More constraints on the 
parameter space and other free parameters will be in- 
cluded in the model. 
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